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Abstract 

We extend the model-independent hkelihood analysis of big bang nucleosynthesis 
(BBN) based on ^He and ^Li to allow for numbers of degrees of freedom which differ 
from the standard model value characterized by Ny = 3. We use the two-dimensional 
likelihood functions to simultaneously constrain the baryon-to-photon ratio and the 
number of light neutrinos. The upper limit thus obtained is N^, < 4.0 (at 95% C.L.). 
We also consider the consequences if recent observations of deuterium in high-redshift 
QSO absorption-line systems (QSOALS) are confirmed. 
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1 Introduction 



In recent years the comparison between standard big bang nucleosynthesis and observational 
data on elemental abundances has come under close scrutiny. Minor changes in the theoret- 
ical abundances (improvements in the code, updated neutron lifetime data) combined with 
many new determinations of the abundances of the light element isotopes, D, ^He, '^He, and 
^Li have led to a flurry of sophisticated statistical analyses of BBN |^, ^, ^, ^ |^, |^ || . The 
results of these analyses however depend strongly on the assumptions made about galactic 
chemical evolution, and in particular, the chemical evolution of deuterium. For the most 
part, the comparison between BBN theory and the observational data made heavy use of 
the solar value of D + ^He because it provided a (relatively large) lower limit for the baryon- 
to-photon ratio 0, rjiQ = lO^^rj > 2.8. This limit for a long time was seen to be essential 
because it provided the only means for bounding rj from below and in effect allows one to 




set an upper limit on the number of neutrino flavors [|T0[, N,y, as well as other constraints 
on particle physics properties. The upper bound to N^, is in fact strongly dependent on 
the lower bound to rj. For lower rj the ^He abundance predicted by BBN drops, allowing 
for a larger N^, to match a given primordial ^He abundance, Yp, as determined from the 
observations. For rj < 4 x 10~^^, there is no bound whatsoever on N^, 

It was argued that since stars (even massive stars) do not destroy ^He in its entirety, 
we can obtain a bound on rj from an upper bound to the solar D and ^He abundances. One 
can in fact limit p, |12| the sum of primordial D and ^He by applying the expression below 
at t = © 

(1) 

In (|l|), is the fraction of a star's initial D and '^He which survives as ^He. For > 0.25 as 
suggested by standard stellar models, and using the solar data on D/H and ^He/H, one finds 
rjiQ > 2.8. This argument has been improved recently ultimately leading to a stronger 
limit ?7io > 3.8 and a best estimate r]io = 6.6 ± 1.4. The stochastic approach used in 
could only lower the bound from 3.8 to about 3.5 when assuming as always that gs > 0.25. 

The limit rjio > 2.8 derived using (|I|) is really a one shot approximation. Namely, it is 
assumed that material passes through a star no more than once. However, to determine 
whether or not the solar (and present) values of D/H and ^He/H can be matched to a given 
primordial abundance, it is necessary to consider models of galactic chemical evolution. In 
the absence of stellar ^He production, particularly by low mass stars, it was shown that 
there are indeed suitable choices for a star formation rate and an initial mass function to: 
1) match the D/H evolution from a primordial value (D/H)p = 7.5 x 10~^, corresponding to 
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7]iQ = 3, through the solar and present interstellar medium (ISM) abundances, while 2) at 
the same time keeping the ^He/H evolution relatively flat so as not to overproduce '^He at 
the solar and present epochs. This was achieved for > 0.3. 

In the context of models of galactic chemical evolution, there is, however, little justifica- 
tion a priori for neglecting the production of ^He in low mass stars. Indeed, stellar models 
predict that considerable amounts of ^He are produced in stars between 1 and 3 Mq [1^ and 



this prediction is consistent with the observation of high ^He/H in planetary nebulae |17]. 
Generally, implementation of these ^He yields in chemical evolution models leads to a gross 
overproduction of ^He/H particularly at the solar epoch |T^ . 

As indicated earlier, the presence of a lower bound on r] allows us to place an upper 
bound to the number of neutrino flavors. From (|l]), the limit rjio > 2.8 corresponds to the 



limit Nu < 3.3 []T9. However, it should be noted that for values of rjio larger than 2.8, the 



central or best-fit value for is closer to 2 [g^, y, |[ and the upper bound is actually found 
to be much smaller with a careful treatment of the uncertainties, < 3.1 [EH, R], though 



this limit is relaxed somewhat when the distribution for is renormalized ||2T|. The range 
in r/io of 6.6 ± 1.4, corresponds to an even tighter limit on = 2.0 ± 0.3 |^ and indicates 
a problem when trying to make use of D and ^He in conjunction with ^He. 

Given the magnitude of the problems concerning ^He |]18|, it would seem unwise to make 
any strong conclusion regarding big bang nucleosynthesis which is based on ^He. In addition, 
deuterium is highly sensitive to models of chemical evolution. Values of primordial D/H as 
high as 2 X 10~^ (i.e. as determined in some measurements of D/H in quasar absorption 
systems [^) were shown to be viable in models of chemical evolution and consistent with 
other observational constraints ||2^. For these reasons it was argued |0, H that the analysis 



comparing BBN theory and observations should be based primarily on the two isotopes 
which are the least sensitive to the effects of evolution, namely, ^He and ^Li. 

In P], the baryon-to-photon ratio, rj was determined on the basis of a likelihood analysis 
using ^He and ''Li. Because of the monotonic dependence of ^He on rj, the ^He likelihood 
distribution has a single (but relatively broad) peak at r^io = 1.75. ''Li, on the other hand, 
is not monotonic in 7], the BBN prediction has a minimum at r/io — 3 and as a result, for an 
observationally determined value of ''Li above the minimum, the ^Li likelihood distribution 
will show two peaks, in this case at r]iQ = 1.80 and r/io = 3.5. The combined likelihood distri- 
bution reflecting both '^He and ''Li is simply the product of the two individual distributions. 
It was found that when restricting the analysis to the standard model, including = 3, 
that the best fit value for tjiq is 1.8 and that 

1.4 < riio < 4.3 95%CL (2) 

In determining (0) systematic errors were treated as Gaussian distributed. When D/H from 
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quasar absorption systems (those showing a high value for D/H |^2[) is included in the 
analysis this range is cut to 1.50 < r^io < 2.55. 

In 1^, the range in t] from the likelihood analysis was translated into a range for N^. 
Because of the agreement between the likelihood predictions of ^He and ^Li for N^, = 3, the 
best fit value was found to be very close to A^,^ = 3. Since the predicted "^He abundance is 
sensitive to N^, (and 77), the uncertainties in in the helium abundance can be translated into 
an uncertainty in N^. The resulting best-fit to N^, based on ^He and ^Li was found to be @] 

K = 3.0 ± 0.2 ± 0.4+°;^ (3) 

thus preferring the standard model result of = 3 and leading to N^, < 3.90 at the 95 % CL 
level when adding the errors in quadrature. In (|^), the first set of errors are the statistical 
uncertainties primarily from the observational determination of Y and the measured error in 
the neutrino half-life, r„. The second set of errors is the systematic uncertainty arising solely 
from ^He, and the last set of errors from the uncertainty in the value of rj and is determined 
by the combined likelihood functions of ^He and ''Li, i.e. taken from Eq. (|^). A similar result 
is obtained when the high D/H value indicated by some observations of quasar absorption 
systems is included in the analysis. 

In this paper, we follow the approach of |^ and ^ in constraining the theory solely 
on the basis of the ^He and ^Li data. We extend the approaches of those two papers by 
fully exploring the {N^,, rf) parameter space of the theory in a self consistent way. We also 
explore the consequences of the QSO deuterium data. In the next section we discuss the 
observational data, and explain our choices of "observed abundances" . Section 3 covers the 
likelihood functions we use. We discuss our results in section 4, and draw conclusions in 
section 5. 



2 Observational Data 

The most useful site for obtaining ^He abundances has proven to be H ll regions in irregular 
galaxies. Such regions have low (and varying) metallicities, and thus are presumably more 
primitive than such regions in our own Galaxy. Because the ^He abundance is known for 
regions of different metallicity one can trace its evolution as a function of metal content, and 
by extrapolating to zero metallicity we can estimate the primordial abundance. 

There is a considerable amount of data on ^He, O/H, and N/H in low metallicity ex- 
tragalactic H ll regions |25|, E6l. In fact, there are over 70 such regions observed with 



metallicities ranging from about 2-30% of solar metallicity. The data for ^He vs. either O/H 
or N/H is certainly well correlated and an extensive analysis pO] has shown this correlation 
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to be consistent with a linear relation. While individual determinations of the '^He mass 
fraction Y have a fairly large uncertainty {AY > 0.010), the large number of observations 
lead to a statistical uncertainty that is in fact quite small PUI E^. Recent calculations [53] 



based on a set of 62 metal poor (less than 20 % solar) extragalactic H ll regions give 

Fp = 0.234 ± 0.002(stat.) ± 0.005(syst.) (4) 
When this set is further restricted to include only the 32 H ll with oxygen abundances less 



than 10% solar it was found that 28 



Yp = 0.230 ± 0.003(stat.) ± 0.005(syst.) (5) 

We will consider both values in our analysis below. 

The primordial ^Li abundance is best determined by studies of the Li content in various 
St function of metallicity (in practice, the Fe abundance). At near solar metallicity, 

the Li abundance in stars (in which the effects of depletion are not manifest) decreases with 
decreasing metallicity, dropping to a level an order of magnitude lower in extremely metal 
poor Population II halo stars with [Fe/H] < —1.3 ([Fe/H] is defined to be the log^^Q of the 
ratio of Fe/H relative to the solar value for Fe/H). At lower values of [Fe/H], the Pop II 
abundance remains constant down to the lowest metallicities measured, (some with [Fe/H] 



< —3 !) and form the so-called "Spite plateau" [|2^. With Li measured for nearly 100 such 
stars, the plateau value is well established. We use the recent results of [^ to obtain the 
^Li abundance in the plateau 

^Li 

1/7 = — = (1.6 ± 0.07) X 10"^° (6) 
H 

where the error is statistical. Again, if we employ the basic chemical evolution conclusion 
that metals increase linearly with time, we may infer this value to be indicative of the 
primordial Li abundance. 

One should be aware that there are considerable systematic uncertainties in the plateau 
abundance. First there are uncertainties that arise even if one assumes that the present Li 
abundance in these stars is a faithful indication of their initial abundance. The actual ''Li 
abundance is dependent on the method of deriving stellar parameters such as temperature 
and surface gravity, and so a systematic error arises due to uncertainties in stellar atmo- 
sphere models needed to determine abundances. While some observers try to estimate these 
uncertainties, this is not uniformly the practice. To include the effect of these systematics, 
we will introduce the asymmetric error range Ai = li^''^ which covers the range of central 
values for ''Li/H, when different methods of data reduction are used (see eg. [p^). Another 
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source of systematic error in the ^Li abundance arises due to uncertainty as to whether the 
Pop II stars actually have preserved all of their Li. While the detection of the more fragile 
isotope ^Li in two of these stars may argue against a strong depletion it is difficult 



to exclude depletion of the order of a factor of two. Furthermore there is the possibility 
that the primordial Li has been supplemented, by the time of the Pop II star's birth, by 
a non-primordial component arising from cosmic ray interactions in the early Galaxy ||32|| . 
While such a contribution cannot dominate, it could be at the level of tens of percent. We 
considered the effects of these corrections in and will not pursue them further here. 

Finally, there have been several recent reported measurements of D/H in high redshift 
quasar absorption systems. Such measurements are in principle capable of determining the 
primordial value for D/H, and hence f], because of the strong and monotonic dependence of 
D/H on 1]. However, at present, detections of D/H using quasar absorption systems indicate 
both a high and low value of D/H. As such, it should be cautioned that these values may 
not turn out to represent the true primordial value. The first of these measurements 
indicated a rather high D/H ratio, D/H ^ {1.9 - 2.5) xlO~^. A re-observation of the high 
D/H absorption system has been resolved into two components, both yielding high values 



with an average value [33 



1/2 = D/H = (L9 ± 0.4) X 10"^ (7) 



Other high D/H ratios were reported in |3J]. However, there are reported low values of D/H 
in other such systems |35] with values D/H ~ 2.5 x 10~^, significantly lower than the ones 
quoted above. 

In addition to our analysis based on ^He and ^Li, we will also consider as in the effects 
of including D/H in the likelihood analysis on rj and N^. We will present results which make 
use of the high D/H measurements. We will comment on the implications of the low D/H 
measurements on our results as well. 



3 Likelihood Functions 

Monte Carlo and likelihood analyses have by now become a common part of the study of 
BBN H, I, I, I, I, I, 0, 1. Our likelihood analysis follows that of [g, except that we now 
consider the additional degree of freedom, A^i,. We begin with a likelihood function for the 
predicted value of Yp 

L,,T^eorAy, N., V) = ^ , exp i^^^fM^^] (8) 

V2^ay(Ar^,r/) V ^(rpiK^T]) J 
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where Yp{Nu,ri) and aY{N^,ri) represent the results of the theoretical calculation. In g], 
we treated the systematic errors in the '^He observations as a either a flat distribution, or 
as a gaussian. This led to L4 0bs(^) being represented by the difference between two error 
functions, or as a gaussian, respectively. Here, we follow the latter prescription, combining 
the statistical and systematic errors in Yp by adding them in quadrature. The resulting form 
for L4 obs(^) is gaussian and can be expressed as 

^4,obs(l^) = exp i^^r^] (9) 

where Yq and ayo characterize the observed distribution and are taken from Eqs. and 
(1^). A full likelihood function for ^He is then obtained by convolving the two likelihood 
functions representing the theory and observational data on ^He 

HNu, V)= JdY L4,Theory(>^, N,, v)U,Ohs{Y) (10) 

This then leads to a complete likelihood function for '^He 

We can construct likelihood functions for ^Li and D in a similar manner, again adding 
statistical and systematic errors in quadrature 

L^iK. V) = ^ exp ( -M^->^) -^:o);) (12) 

^ " v/27r(aKiV.,r/)+a?o) ^ l2(^?(^- ^) + ^?o) ^ ^ ' 

UK. V) = ' exp ( -M^->^) -l/2o);^ . (13) 

V'2vr(ai(iV.,r^)+aio) ' \ 2{aUK,v) + ^^lo) J ^ ' 

In these expressions, y^Q and ?/2o are the observed values of the ^Li and D abundances, and a-jQ 
and (720 their associated uncertainties, while those quantities which are shown as functions 
of {Nu^T]) are all derived from our theoretical BBN calculations. 

The quantities of interest in constraining the Ny-rj plane are the combined likelihood 
functions 

L47 = L4 X L7 (14) 

and 

L247 = L2 X L47. (15) 
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Contours of constant L47 represent equally likely points in the N,y-ri plane. Calculating 
the contour containing 95% of the volume under the L47 surface gives us the 95% likelihood 
region. If we wish to impose a constraint in the analysis (for example, > 3) we can search 
for the contour which contains 95% of the volume under L47 which satisfies the constraint. 
From these contours we can then read off ranges of N^, and 77. 

4 Results 

Using the abundances in eqs (^J^, and |^ and adding the systematic errors to the statistical 
errors in quadrature we have 

Fp = 0.234 ± 0.0054 (16) 

1/7 = (1.6 ±0.36) X 10^1° (17) 

1/2 = (1.9 ±0.4) X 10-^ (18) 

These abundances are then combined with our theoretical calculations to produce the like- 
lihood distributions discussed above. 

Figure 1 shows a three-dimensional view of both L47 and L247 (in arbitrary units). As 
one can see, L47 is double peaked. This is due to the minimum in the predicted lithium 
abundance as a function of t] as was discussed earlier. The peaks of the distribution as well 
as the allowed ranges of 77 and N^, are more easily discerned in the contour plot of Figure 2 
which shows the 50%, 68% and 95% confidence level contours in L47 and L247. The crosses 
show the location of the peaks of the likelihood functions. Note that L47 peaks at N,^ = 3.0, 
Tjio = 1.8 (in agreement with our previous results §]) and at N^, = 2.3, rjio = 3.6. The 95% 
confidence level allows the following ranges in 77 and N^, 

1.6 <N^< 4.0 

1.3 < 7/10 < 5.0 (19) 

Note however that the ranges in rj and iVj, are strongly correlated as is evident in Figure 2. 

Since L2 picks out a small range of values of r/, largely independent of N^, its effect on 
L247 is to eliminate one of the two peaks in L47. L247 also peaks at N^, = 3.0, r/10 = 1.8. In 
this case the 95% contour gives the ranges 

2.0 <N^< 4.1 

1.4 < 7/10 < 2.6 (20) 
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Note that the additional constraint has raised the upper hmit on (shghtly) though this 
is counter to intuition. 

These results can be compared with those of who have used a Bayesian analysis to 



examine BBN. Though they have not performed the full two-dimensional likelihood analysis 
presented here, their analog of L47 peaks at N^, — Ay/0.016 ~ 2.8 and 3.9. Since they have 
taken a significantly larger value for the observed ^He abundance (Yp = 0.242 + AF ib 0.003) 
this corresponds to Ny k, 2.3 and 3.4 (using AY = —0.008). The remaining difference 
between our results is presumably due to the difference in our assumptions about the observed 
^Li abundance. Our results for L247 also appear to be consistent with those of ^7 . 



Since Njy = 3 is well within the range covered by the 95% contour, it is legitimate to 



impose the additional constraint A^^ > 3 Figure 3 shows the 50%, 68% and 95% 



confidence level contours under this condition. The unconstrained contours are also shown 
(dashed) for comparison. This has a rather small effect on the upper limit to N^. In the 
case of L47 it has increased to N,, < 4.1, while for L247 it has dropped to N^, < 4.0. 

Finally, in figure 4 we consider our alternative choice for the ^He abundance based on 
the lowest metallicity set of H ll regions [^] from eq. (^ which when adding the systematic 



errors in quadrature to the statistical errors gives 

= 0.230 ± 0.0058. (21) 

The solid curves show the confidence level contours under these circumstances (the dashed 
curves show the old contours, for comparison). As expected, the lower value of Yp drives all 
the contours towards lower values of N^. The peaks occur at N^, = 2.7, rjiQ = 1.7 and at 
Nu = 2.1, rjio = 3.4. The corresponding ranges are, for L47 (unconstrained): 1.3 < Ny < 3.7, 
1-3 < TjiQ < 5.0; for L247 (unconstrained): 1.8 < A^^ < 3.9, 1.3 < r^io < 2.5. When the 
constraint > 3 is applied, we find for L47 {Ny > 3): A^^ < 4.0 and for L247 (A^^ > 3): 
Ny < 3.9. 



5 Conclusions 

We have presented a full two-dimensional likelihood analysis based on big bang nucleosyn- 
thesis and the observations of ^He,''Li and certain (high) determinations of D/H in quasar 
absorption systems. Allowing for full freedom in both the baryon-to-photon ratio, rj, and the 
number of light particle degrees of freedom as characterized by the number of light, stable 
neutrinos, A^^,, we have confirmed the successful predictions of BBN in a limited range in 77 
and a range in Ny which not only encompasses the standard model value of Ny = 3, but 
whose likelihood is peaked at or near that value. The largest value of Ny allowed at the 95% 
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CL is found to be 4.0 at a value of t]iq = 1.8, when using only ^He and ''Li in the analysis. 



Because the high values of D/H observed in certain quasar absorption systems [22, |3^, Q 
agree so well with the predictions made when using ^He and ''Li this result is only slightly 
altered when including D/H. In contrast, had we used the low D/H values found in ||35|| , 
and corresponding to a central value of rjiQ ~ 6.4, there would be virtually no overlap in the 
likelihood distributions of this value of D/H and L47 [^. Though we would not claim that 



these results indicate a problem with the low D/H measurements per se, they do indicate 
the degree to which they are incompatible with the current determinations of ^He and ^Li. 
Until this incompatibility is resolved, no useful limit on can be derived with this method 
when the low D/H abundance is used. 
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Figure 1: Likelihood functions L47 (top) and L247 (bottom) using Yp = 0.234 ± 0.0054. 
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Figure 2: The top panel shows contours in the hkehhood function L47 for Yp = 0.234±0.0054. 
The contours represent 50% (innermost), 68% and 95% (outermost) confidence levels. The 
crosses mark the points of maximum likelihood. The bottom panel shows the equivalent 
results for L247 (with L47 shown dashed, for comparison.) 
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Figure 3: The top panel shows the hkehhood function Z/47 given the constraint N^, > 3. The 
contours represent 50%, 68% and 95% confidence levels. The unconstrained likelihood is 
shown dashed. In both cases Yp = 0.234 ± 0.0054. The bottom panel shows the equivalent 
results for I/247. 
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Figure 4: The top panel shows hkehhood contours for L47, with Yp = 0.230 ± 0.0058. The 
contours represent 50%, 68% and 95% confidence levels. Results for Yp = 0.234 ±0.0045 are 
shown dashed for comparison. The bottom panel shows the equivalent results for L247. 
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